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ABSTRACT 

The  problem  of  simultaneous  determination  of  optimal  due  dates  and  optimal  sequence  for 
N-job  single-machine  problem  with  multiple  due  dates  is  considered  in  this  paper.  The  penalty  for 
a  job  is  assumed  to  be  a  linear  function  of  the  due  date  and  the  earliness/tardiness  for  the  job.  The 
objective  is  to  minimize  the  total  penalty  for  all  jobs.  An  efficienf  optimal  algorithm  to  solve  the 
problem  is  developed  and  several  important  results  are  proved. 
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INTRODUCTION  •  ..     - 

This  paper  considers  the  problem  of  simultaneous  determination  of  optimal  due  dates  and 
optimal  sequence  for  N  independent  jobs  to  be  processed  on  a  single  machine.  All  jobs  are 
available  at  time  zero.  Processing  times  are  known  and  deterministic.  The  penalty  for  a  job  in  a 
sequence  is  assumed  to  be  a  linear  function  of  the  due  date  assigned  to  the  job  and  the 
earliness/tardiness  for  the  job  in  the  sequence.  The  same  linear  penalty  function  is  used  for  all 
jobs.  The  objective  is  to  minimize  the  total  penalty  for  all  jobs.  The  number  of  distinct  due  dates, 
m,  to  be  assigned  to  the  jobs  is  assumed  to  be  prespecified  and  known.  We  allow  m  to  take  any 
value  in  the  set  {1,2,...,N}.  Job  preemption  is  not  allowed. 

The  problem  of  simultaneous  determination  of  optimal  due  dates  and  optimal  sequence  has 
been  motivated  and  considered  by  several  authors  in  the  recent  literature  [1,  2,  6,  9,  10];  the  reader 
is  directed  to  Baker  and  Scudder  [3]  and  Bector,  Gupta  and  Gupta  [4]  for  a  brief  review  of  this 
literauire.  The  problem  considered  in  this  paper  is  most  directly  related  with  the  problems 
considered  in  Panwalkar,  et  al.  [5],  Seidman,  et  al.  [7],  and  Bagchi  [1].  Panwalkar,  et  al.  [5] 
assume  a  common  due  date  for  all  jobs;  that  is,  m  =  1 .  The  penalty  function  used  in  this  paper  is 
the  same  as  in  [5].  Several  results  from  [5]  are  used  in  this  paper.  As  in  [5],  we  show  that  the 
penalty  function  can  be  expressed  as  a  product  of  two  sequences,  which  is  minimized  by  arranging 
them  in  opposite  order.  Seidman,  et  al.  [7]  consider  the  problem  when  each  job  is  allowed  a 
different  due  date.  They  show  that  an  SPT  (shortest  processing  time)  sequence  is  optimal  for  this 
problem.  With  the  number  of  distinct  due  dates,  m,  permitted  to  be  any  number  in  the  set 
{ 1,2,. ..,N},  our  paper,  in  essence,  bridges  the  gap  between  the  scheduling  models  with  a 
common  due  date  and  the  scheduling  models  with  possibly  N  distinct  due  dates  for  the  N-job 
problem. 

Bagchi  [1]  assumes  that  the  N  jobs  to  be  processed  on  the  single  machine  are  to  meet  the 
requirements  of  m  customer  orders.  It  is  assumed  that  m  <  N  and  that  each  customer  order  has  one 
or  more  jobs  in  it.  The  specific  jobs  in  each  customer  order  are  assumed  to  be  known.  A  common 
due  date  is  assumed  for  all  jobs  in  a  customer  order.  Thus,  there  could  be  m  distinct  due  dates  for 


all  jobs.  Instead  of  the  due  date  penalty  used  in  our  paper,  Bagchi  uses  a  lead  time  penalty  which 
is  a  linear  function  of  the  lead  time,  where  lead  time  for  a  job  in  a  customer  order  is  the  completion 
time  of  the  last  job  in  the  customer  order.  TTiese  assumptions  allowed  Bagchi  to  decompose  the 
sequencing  problem  into  two  independent  sequencing  subproblems  -  one  to  schedule  the  jobs 
within  a  customer  order  and  the  other  one  to  schedule  the  customer  orders.  The  first  subproblem  is 
shown  to  be  the  common  due  date  problem  of  Panwalkar,  et  al.  [5].  The  second  subproblem  is 
shown  to  be  the  minimum  weighted  completion  time  problem  which  can  be  solved  using  the 
WSPT  (weighted  shortest  processing  time)  algorithm. 

Our  paper  considers  two  different  cases.  In  the  first  case  we  assume  that  the  decision 
maker  specifies  the  number  of  jobs  to  be  assigned  to  different  due  dates.  However,  unlike  Bagchi 
[1],  the  specific  jobs  to  be  assigned  to  different  due  dates  are  unknown.  In  the  second  case,  we 
assume  that  the  number  of  jobs  to  be  assigned  to  different  due  dates  is  unknown.  It  should  be 
remarked  that  the  second  case  turned  out  to  be  very  difficult  to  analyze. 

Section  2  of  the  paper  defines  the  notation  and  formulates  the  problem. 

Section  3  of  the  paper  considers  the  case  when  the  vector  (ni,n2,...,nm),  where  nj  is  the 
number  of  jobs  assigned  to  the  jth  due  date,  is  externally  specified.  Note  that  while  nj  is  assumed 
to  be  known,  the  specific  jobs  to  be  assigned  to  the  jth  due  date  are  unknown  and  to  be  determined. 
It  is  shown  that  the  objective  function  can  be  expressed  as  a  sum  of  the  pairwise  product  of  two 
equal  sets  of  numbers,  where  one  set  consists  of  processing  times  and  the  other  set  consists  of 
positional  weights.  The  positional  weights  are  independent  of  job  characteristics.  An  efficient 
procedure  is  presented  to  find  optimal  due  dates  and  an  optimal  sequence.  The  results  in  this 
section  are  used  in  the  next  section  in  developing  an  optimal  procedure  when  the  vector 
(ni,n2,...,nm)  is  also  a  decision  variable. 

Section  4  of  the  paper  develops  some  dominance  results  to  find  an  optimal  vector 
(ni,n2,...,nm)  when  it  is  not  externally  specified.  These  dominance  results  are  used  in  developing 
an  optimal  procedure  to  find  (ni,n2,...,nm).  Even  though  these  results  are  intuitive,  proofs  for 
these  turned  out  to  be  very  difficult. 


Section  5  of  the  paper  discusses  the  computational  requirements  for  the  algorithm  in 
Section  4.  It  is  shown  that  our  algorithm  can  give  large  computational  savings  compared  to  a 
complete  search.  In  one  example,  we  show  that  our  algorithm  is  more  than  lO^  times  faster  than  a 
complete  search. 

The  paper  closes  in  Section  6  by  providing  some  concluding  remarks  and  ideas  for  further 
research. 

2.  PROBLEM  FORMULATION 

Throughout  the  paper,  we  assume  that  the  number  of  jobs  to  be  scheduled  and  the  number 
of  due  dates  to  be  assigned  to  the  jobs  are  known.  The  following  notation  is  used  to  formulate  the 
problems: 

N      =      Number  of  jobs  to  be  scheduled, 

m      =      Number  of  due  dates, 
(n,n+k>       =      {n,n+l,n+2,...,n+k}, 

di      =      Due  date  for  job  i,  i  G  <1,N>, 

Q       =      Completion  time  for  job  i,  i  G  (1,N>, 

Ei       =      Max  (0,di-Ci)  =  Earliness  for  job  i,  i  G  (1,N>, 

Ti       =      Max  (0,Ci-di)  =  Tardiness  for  job  i,  i  g  <1,N>, 
ti       =      Processing  time  for  job  i,  i  g  (1,N>, 

Pi       =      Per  unit  time  due  date  penalty  for  each  job, 

P2       =      Per  unit  time  earliness  penalty  for  each  job,  and 

P3       =      Per  unit  time  tardiness  penalty  for  each  job. 
In  addition,  we  let  D]  <  D2  <  ...  <  D^  denote  the  m  due  dates  and  let  Ij  denote  the  set  of  jobs 
assigned  to  due  date  Dj  for  j  g  (l,m).  As  in  Panwalkar,  Smith  and  Seidman  [5],  the  constants  Pi, 
P2  and  P3  are  assumed  to  be  non-negative  and  known. 

The  objective  is  to  determine  D  =  {Di,D2,...,Dm},  I  =  {Ii,l2,...,Im}  and  a  schedule  a  to 
minimize  the  total  penalty  TP(T),I,a),  where 


TP(D,I,o)=   I     I    {Pi-Dj+P2-Ei+P3-T,}  (1) 

J=MeIj 
and  di  =  Dj  for  i  g  Ij. 

It  is  easy  to  see  that  the  problem  has  an  optimal  solution  with  zero  machine-idle  time.  (See 

Lemma  2  in  the  Appendix.)  Therefore,  it  is  sufficient  to  consider  permutation  schedules  to  fmd 

optimal  solutions. 

Let  [k]  denote  the  job  in  position  k  in  a  permutation  schedule,  then 
k 

C[k]=It[i].  (2) 

1=1 

We  will  use  a  to  denote  a  permutation  schedule  in  the  rest  of  the  paper.  Also,  jobs  are  assumed  to 

be  indexed  such  that 

0<ti<t2<...<tN  (3) 

The  next  section  considers  the  problem  when  lljl  =  nj  is  known  for  all  j  e  (l,m>.  Note  that 

while  lljl  is  assumed  to  be  known,  the  specific  jobs  in  Ij  are  unknown  and  to  be  determined. 


3.  RESULTS  AND  PROCEDURE  WITH  lljl  =  nj  KNOWN 

This  section  assumes  that  the  vector  (ni,n2,...,nm)  of  m  positive  integers  such  that 

m 

X  nj  =  N  and  lljl  =  nj  for  j  =  1,2,. ..,m  is  externally  specified.  D,  I  and  a  are  decision  variables. 
j=l  '  ^        ' 

The  cases  when  m=l  and  m=N  are  considered  in  Panwalkar,  Smith  and  Seidman  [5]  and  Seidman, 

Panwalkar  and  Smith  [7],  respectively. 

For  Pi  >  P3,  an  SPT  sequence  is  optimal  with  dj  =  0  for  i  e  (1,N>  (see  Theorems  1  and  2 

in  [7]).  We  assume  Pi  <  P3  in  developing  the  properties  and  the  algorithm  below. 


Properties  of  Optimal  Solutions:  We  now  state  several  properties  of  optimal  solutions  which  will 

be  used  in  developing  an  optimal  algorithm.  Lemmas  1  and  3  needed  to  develop  these  properties 

j 
are  proved  in  the  appendix.  We  define  Nj  =   X  "k  for  j  e  <l,m>  with  Nq  =  0;  Nj  is  the  total 

k=l 

number  of  jobs  assigned  to  the  first  j  due  dates. 


From  Lemma  1  in  the  Appendix,  it  is  easy  to  see  that  for  any  given  D  and  a,  there  is  an 
optimal  I  such  that  a  batch  of  consecutive  nj  jobs  in  o  is  assigned  to  Dj.  The  following  property 
gives  an  optimal  I  for  any  given  D  and  a. 

Propertv  1:  For  any  given  D  and  o,  there  is  an  optimal  I  such  that 

^j  =  {  ^Nj.i  +  1'  ^Nj.i+2'  •••'  ^Nj.i+nj  j 
for  j  e  (l,m),  where  an  is  the  job  in  position  n  in  sequence  a. 

This  property  essentially  says  that  there  is  an  optimal  solution  such  that  nj  consecutive  jobs 
(in  positions  Nj.i+1  to  Nj.i+nj)  in  a  are  assigned  to  Dj. 

Lemma  3  in  the  Appendix  shows  that  for  any  given  a,  there  is  an  optimal  D  such  that  Dj 
coincides  with  the  completion  time  of  a  job  in  o.  The  following  property,  which  gives  an  optimal 
D  for  a  given  a,  follows  from  Panwalkar,  et  al.  [5]. 

Property  2:  For  a  given  a,  there  is  an  optimal  D  such  that  Dj  =  Cfk;],  where 

kj  =  Nj.,+(^^-njj*  (4) 

(x)+  here  stands  for  the  smallest  integer  >  x. 

Property  2  implies  that  the  first  (kj-Nj.i)  jobs  in  Ij  are  early  and  the  remaining  (Nj-kj)  jobs 

are  tardy. 

We  now  develop  an  algorithm  to  find  (D,I,a)  to  minimize  (1). 

N 
Algorithm:  Using  Properties  1  and  2,  we  first  reduce  (1)  to  the  form  £  wi'tfii.   Knowing  the 

i=l 
positional  weights  wj's  for  positions  i  e  (1,N),  an  optimal  sequence  can  be  found  as  in 

Panwalkar,  Smith  and  Seidman  [5]. 

Using  Property  1,  we  get: 

m         Nj 
TP(D,a)  =  Min  TP(D,o,I)  =   S  t  (PrDj+P2-E[i]+P3-T[i]). 

I  j=l    i=Nj.i  +  l 

Further,  using  Property  2,  we  can  write: 


m 
TP(a)  =  Min  TP(D,o)  =  X 
D  j=l 


Ni 
nj.PrCi 


:[k.]+.    1  P2-(C[k  ]-C[i])+.  rP3-(C[i]-C[kj]) ,  (5) 

■'     i=Nj.i  +  l        ••  i=kj+l  ■'   . 


k 
with  kj  defined  by  (4).  Using  C[k]  =  Z  t[i]  for  a  given  a,  (5)  can  be  simplified  to: 
N  i=l 

TP(a)=  Iwit[i],  (6) 

i=l 
where 

fP2-(i-l-Nj.i)+Pr(N-Nj.i)        for  i  g  <Nj.i  +  l,kj> 
,      Wi  =  <  (7) 

[P3-(Nj-i+l)+Pi-(N-Nj)  foriG  <kj+l,Nj> 

Note  that  (6)  is  a  product  of  two  sequences,  which  is  minimized  by  arranging  them  in 

opposite  order.  The  following  steps  give  an  optimal  solution  for  a  given  (ni,n2,.",nni).  We  call 

this  algorithm  the  Given-(ni,n2,.",nm)-Algorithm  or  the  Gn-Algorithm. 


Steps  of  the  Gn-Algorithm 

J 
Step  0:     Set  No  =  0  and  Nj  =   I  nk  for  j  e  <l,m>. 

k=l 
Step  1:     For  j  s  (l,m),  compute 

kj  =  Nj.,+(^.nj)'. 
Step  2:     For  i  g  (1,N>,  compute  w,  using  (7). 
Step  3:     Rank  the  positional  weights  wi's  in  a  descending  order  of  magnitude  such  that  the 

largest  w,  is  ranked  first  and  the  smallest  wi  is  ranked  N^.  Break  the  des 

arbitrarily. 
Step  4:     Find  a  sequence  o  by  assigning  job  i  to  position  k  where  k  is  such  that  wk  has 

rank  i.  Set 

Dj  =  I  t[i],  and 

1=1 
di  =  Dj  for  i  G  Ij  =  [  o^.^^^ ,  aN^^^2'  -'  ^Nj  } 

for  j  e  (l,m>.  Stop,  the  solution  found  in  Step  4  is  an  optimal  solution. 
It  is  easy  to  see  that  the  first  (kj-Nj.i)  jobs  in  Ij  are  in  an  LPT  order  while  the  last  (Nj-kj) 
jobs  are  in  an  SPT  order.  We  now  give  a  numerical  example  to  illustrate  the  procedure. 


Numerical  Example 

Given  10  jobs  with  ti  =  3,  t2  =  10, 13  =  1 1, 14  =  13, 15  =  13,  t6  =  16, 17  =  17,  tg  =  20, 
t9  =  22  and  tio  =  25.  The  penalties  are  Pi  =  2,  P2  =  1 1  and  P3  =  18.  We  are  given  m  =  2,  ni  =  4 
and  n2  =  6. 

From  Steps  0  and  1  of  the  algorithm,  we  get  Ni  =  4,  N2  =  10,  ki  =  3  and  k2  =  8.  The  ten 
positional  weights  and  their  ranks  are 


Position 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Weight  wi 

20 

31 

42 

30 

12 

23 

34 

45 

36 

18 

Rank 

8 

5 

2 

6 

10 

7 

4 

1 

3 

9 

Optimal  Sequence:  8-5-2-6-10-7-4-1-3-9;  Di  =  t8+t5+t2  =  43,  and  D2  = 
t8+t5-(-t2+t6+tio+t7+U+ti  =  117;  Ii  =  {8,5,2,6}  and  I2  =  {10,7,4,1,3,9};  Penalty  =  3763. 
Next  section  considers  the  problem  when  the  vector  (ni,n2,...,nni)  is  unknown. 

4.  RESULTS  AND  PROCEDURE  WITH  lljl  =  nj  UNKNOWN 

This  section  considers  the  problem  when  the  vector  (ni,n2,...,nm)  is  unknown.  D,  I,  and 
G  are  decision  variables  in  this  case.  Dominance  properties  are  developed  to  find  optimal 
ni,n2,...,nm.  A  key  result  in  the  section  shows  that  there  is  an  optimal  solution  with 
ni  >  n2  >  ...  >  nm.  This  decreasing  property  is  used  to  develop  an  efficient  search  algorithm  to 
find  optimal  ni,n2,...,nni.  Computational  requirements  for  this  search  algorithm  are  discussed. 

Recall  that  w,  denotes  the  positional  weight  associated  with  position  i  for  a  known  vector 
(ni,n2,...,nm).  Step  3  of  the  Gn-Algorithm  gives  ranks  for  these  weights.  We  need  the  following 
additional  notation  associated  with  these  weights: 

tti  =  the  positional  weight  with  the  rank  of  i 

N 
WTi=  Sttk 

k=i 

Clearly  ai  >  a2  >  ...  >  a^. 

Note  that  the  weights  ai's  and  WTi's,  while  they  depend  on  the  vector  (ni,n2,..Mnm)  and  the 

penalties  Pi,  P2  and  P3,  they  do  not  depend  on  the  processing  times  for  the  jobs. 


8 

To  develop  the  results  in  this  section,  we  first  assume  m  =  2.  We  let  WTi(n)  denote  the 

weight  WTi  when  ni  =  n  and  n2  =  N-n.  ai(n)  is  similarly  defined.  Let  COST(n)  denote  the 

corresponding  optimum  objective  function  value,  then  we  have 

COST(n)     =      ti-ai(n)  +  t2-a2(n)  +  ...  +  tN-aN(n) 

=      tr(WTi(n)-WT2(n))  +  t2(WT2(n)-WT3(n))  +  ...  +  tN-(WTN(n)) 

=      ti-WTi(n)  +  (t2-ti)-WT2(n)  +  (t3-t2)-WT3(n)  +  ...  +  (tN-tN-l)-WTN(n) 

N 
=      ti-WTi(n)+  lAi-WTi(n),  (8) 

i=2 
where  Ai  =  ti-ti-i. 

Lemma  4:  For  m  =  2,  n  g  <1,N-1)  and  P  e  <1,N-1>,  we  have  COST(n)  <  COST(J?)  for  every 
vector  (ti,t2,..-,tN)  satisfying  (3)  if  and  only  if 

WTi(n)<WTi(i?)      forallis  <1,N>  '    '  (9) 


Proof:  From  (8)  we  have 

N 
COST(J?)  -  COST(n)  =  tr(WTi(P)-WTi(n))  +  I  Ai-(WTi(J?)-WTi(n))  (10) 

i=2 
If  (9)  holds  for  given  n  and  Q,  then  COST(P)-COST(n)  >  0  from  (10)  for  any  given  (ti,t2,...,tN) 

satisfying  (3).  If  (9)  does  not  hold  for  all  i  e  (1,N),  then  there  is  a  k  g  <1,N)  such  that 

WTk(n)  >  WTk(i?)  or  WTk(i?)-WTk(n)  <  0.  Then,  we  can  select  a  vector  (ti,t2,...,tN)  satisfying 

(3)  with  Ak  =  tk-tk-i  high  enough  to  make  COST(P)-COST(n)  <  0  in  (10).  This  completes  the 

proof. 

Lemma  5:  For  m  =  2  and  n  g  (oIy]  ^'  ^^  ^^^^  WTi(n)  >  WTil  ry"!   |  for  all  i  e  <1,N>.  (Recall 
that  (x)+  is  the  smallest  integer  >  x.) 


Proof:  A  proof  for  this  is  given  in  the  appendix. 


THEOREM  1:  There  is  an  optimal  solution  (D*,I*,o*)  such  that  n   >  n    >  ...  >  n 


m 


Proof:  First  consider  the  case  m  =  2.  Un\<iy\  ,  then  WTi(n  i )  >  WYA  l-sy-]     for  all  i  e  ( 1  ,N) 
from  Lemma  5,  and  COST(ni)  >  COST!   -j-      from  Lemma  4.  Thus,  there  is  an  optimal  solution 


withn,  >|r^|   > — J — ;  or  n    >  n-. 


The  proof  for  m  =  2  can  be  extended  for  the  m  >  2  case  by  considering  two  adjacent 
batches.  Thus,  there  is  an  optimal  solution  with  n.  >  n.      for  j  =  l,2,...,m-l,  or 


*       *  * 

n   >n    >  ...  >  n 
12  m 


The  following  search  algorithm  gives  an  optimal  (n  ,n  ,. .  .,n   ).  We  call  this  algorithm  the 


r  2' 


m 


Efficient  Search  Algorithm  (or  the  ES-Algorithm). 


Efficient  Search  Algorithm 


r 


For  nk  = 


k-1         N+      ^ 


N-Inj 

-1=1 
m-(k-l) 


k-1     ^ 
N-Snj 

Im-(k-l); 


+  l,...,N-(m-k) 


m-1 
For  nm  =  N-  X  ij 

Find  minimum  total  penalty  for  (ni,n2,...,nm)  using  the  Gn- Algorithm 


Find  (n  ,n-,...,n    )  corresponding  to  the  lowest  total  penalty. 


V  2 


m 


The  following  example  illustrates  the  procedure. 


Numerical  Example:  Consider  the  10-job  problem  with  processing-time  and  penalty  data  used  in 
the  numerical  example  in  Section  3.  We  solve  it  for  m  =  2  and  m  =  3. 
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Solution  for  m  =  2:  To  find  an  optimal  (n    n.),  we  need  to  find  the  minimum  cost  using  the 
Gn- Algorithm  for  ni  =  5,  6,  7,  8  and  9.  Note  that  n2  =  10  -  ni.  The  following  table  gives  the 
costs  using  the  Gn-Algorithm: 


ni 

5 

6 

7 

8 

9 

n-> 

5 

4 

3 

2 

1 

COST 

3,624 

3,586 

3,729 

4,121 

4,738 

We  get  n   =  6  and  n   =  4  corresponding  to  the  minimum  cost  of  $3,586. 


Solution  for  m  =  3:  The  costs  computed  using  the  Gn-Algorithm  for  various  values  of  ni,  n2  and 
n3  in  the  ES-Algorithm  are  given  below. 

ni m [13 COST 

4  3  3  2,757 

4  4  2  2,852 

5  3  2  2,845 

5  4  1  3,023 

6  2  2  3,021 

6  3  1  3,068 

7  2                           1  3,367 
_8 1^ 1 3,829 

We  get  n    =  4,  n.  =  3  and  n    =  3. 

The  next  section  discusses  the  computational  requirements  for  the  ES-Algorithm. 


5.  COMPUTATIONAL  REOUIREMENTS  FOR  THE  EPTTCIENT  SEARCH  ALGORITHM 

The  ES-Algorithm  generates  vectors  (ni,n2,...,nm)  by  partitioning  N  into  m  integer 
m 

numbers  such  that  S  n,  =  N,  and  ni  >  n2  >  ...  >  nm  >  1.  Let  P    denote  the  total  number  of  these 
i=i  1        ^  Ml  j^ 

vectors,  then  the  following  theorem  gives  a  recursive  procedure  to  compute  P^ . 


THEOREM  2:  For  a  positive  integer  n  >  2  and  re  <l,n-l>,  we  have 

p"  =  p"  =  p|  =  i,and 
1         n         I 

p"  =  pn-r  +  pn-r  +    ^  pn-r  (l^^ 

r         1  2  min(r,n-r) 
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Proof:  P"  =  1  because  ni  =  n  is  the  only  feasible  partition.  P"  =  1  because  ni  =  n2  =  ...  =  nn  =  1 
is  the  only  feasible  partition.  We  now  prove  (11). 

We  are  given  r  due  dates  and  n  >  r  jobs  to  be  assigned  to  these  due  dates.  Since  each  due 
date  is  assigned  at  least  one  job,  we  can  first  assign  one  job  to  each  of  the  r  due  dates.  We  are  left 
with  n-r  jobs. 

P"""^  in  the  r.h.s.  of  (1 1)  gives  the  total  number  of  ways  if  the  remaining  n-r  jobs  are 
assigned  to  the  first  due  date.  P^""^  in  the  r.h.s.  of  (1 1)  gives  the  total  number  of  ways  if  the 
remaining  n-r  jobs  are  assigned  to  the  first  two  due  dates  such  that,  out  of  n-r  jobs,  each  due  date 
gets  at  least  one  job  and  the  first  due  date  gets  at  least  as  many  jobs  as  the  second  one.  The 
remaining  terms  in  the  r.h.s.  of  (1 1)  can  be  explained  similarly. 

Computational  Savings:  The  ES-Algorithm  requires  P^  applications  of  the  Gn-Algorithm  to  find 
(n  ,n-,...,n    ).   ( j^.j  )  gives  the  total  number  of  ways  of  assigning  N  jobs  to  m  due  dates  such 
that  each  due  date  gets  at  least  one  job,  so  a  complete  enumeration  to  find  (n  ,n-,...,n    )  would 
require  (  ^.i  )  applications  of  the  Gn-Algorithm.  The  ratio  (  p^.|  )/r    can  be  used  as  a  measure 
of  computational  savings  realized  by  the  use  of  our  dominance  results  in  the  ES-Algorithm.  Table 
1  gives  the  values  of  P^,  (  ^.i )  and  ( j^.j  )/P    for  N  =  40  and  different  values  of  m. 

Results  in  Table  1  indicate  that  large  computational  savings  can  be  realized  by  using  our 
algorithm  compared  to  a  complete  search.  For  N  =  40  and  m  =  20,  Table  1  indicates  that  the 
ES-Algorithm  is  more  than  10^  times  faster  than  a  complete  search.  (Shanthikumar  and  Yao  [8] 
also  noticed  that  the  use  of  the  decreasing  property  narrows  down  the  search  for  an  optimal 
solution  to  a  very  small  subset  of  the  solution  space.) 

Analysis  for  m  =  1  and  m  =  n: 

For  m  =  1,  we  have  n  =  N,  and  only  one  (Py  =  1)  application  of  the  Gn-Algorithm  is 
required.  It  is  easy  to  see  that  our  algorithm  reduces  to  the  algorithm  of  Panwalkar,  Smith  and 
Seidmann  [5]  for  the  m  =  1  case. 
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^  4(  %  XT 

For  m  =  N,  we  have  n^  =  n^  =  ...  =  n    =  1.  For  this  case  also,  only  one  (PT;  =  1) 
application  of  the  Gn-Algorithm  is  required.  It  is  easy  to  see  that 

Nj    =    n*  +  n*  +  ...  +  n*     =    j  for  j  =  1,2,...,N, 

kj     =    J-l+flj^)         =    J  forj  =  1,2,...,N,  and 

wi    =     Pi-(N+l-i)  for  i  =  1,2,.. .,N. 


TABLE  1:  Table  of  Computational  Savings  for  the  ES-Algorithm  for  N  =  40 


m 

pN 

m 

C.) 

c.x 

1 

1 

1 

1 

2 

20 

3f 

1.95 

3 

133 

741 

5.18 

4 

478 

9139 

19.11 

5 

1115 

82,251 

73.76 

6 

1945 

575,757 

296.01 

10 

3590 

2.12  X  10^  (approx.) 

5.90  X  104 

15 

1861 

1.39  X  10 10  (approx.) 

7.46  X  106 

20 

627 

6.36  X  IQlO  (approx.) 

1.01  X  108 

Since  wi  >  W2  >  ...  >  wn,  the  objective  function  is  minimized  by  arranging  the  jobs  in  an 
increasing  order  of  processing  times.  Thus,  an  SPT  sequence  is  optimal  for  this  case. 


A  Negative  Result:  Let  COST(n)  denote  the  minimum  cost  of  the  N-job,  2-due  date  problem  with 
n  jobs  assigned  to  the  first  due  date.  In  the  numerical  example  in  Section  4,  COST(n)  has  the 
incline  property;  that  is,  as  n  increases,  once  COST(n)  goes  up,  it  does  not  come  down.  We  found 
that  COST(n)  has  this  property  for  most  of  the  problems  that  we  solved. 

Unfortunately,  COST(n)  does  not  always  have  this  property  as  shown  by  the  following 
counter-example.  Further  computational  savings  could  be  realized  if  COST(n)  had  the  incline 
property. 
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Counter-Example:  Assume  N  =  6;  Pi  =  17,  P2  =  1 1,  P3  =  18;  ti  =  t2  =  13  =  14  =  t5  =  5,  and  t^ 
25.  The  following  table  gives  COST(n)  for  different  values  of  n. 


n 

1 

2 

3 

4 

5 

COST(n) 

2,195 

2,180 

2,175 

2,180 

2,175 

COST(n)  in  this  example  does  not  follow  the  incline  property. 

6.  CONCLUDING  REMARKS  '■  ' 

This  paper  makes  two  major  contributions.  First,  an  optimal  algorithm  of  the  order 
0(N  log  N)  was  developed  to  determine  optimal  due  dates  and  a  sequence  for  an  N-job  problem 
when  the  number  of  jobs  to  be  assigned  to  different  due  dates  are  known.  It  was  shown  that  the 
objective  function  can  be  expressed  as  a  sum  of  the  pairwise  product  of  two  equal  sets  of  numbers, 
where  one  set  consists  of  processing  times  and  the  other  set  consists  of  positional  weights. 

Second,  an  important  dominance  property  was  developed  to  determine  the  number  of  jobs 
to  be  assigned  to  different  due  dates.  While  the  result  is  intuitive,  the  proof  turned  out  to  be 
difficult.  An  algorithm  of  the  order  0(P^*N  log  N)  was  developed  to  find  optimal  due  dates, 
sequence  and  number  of  jobs  to  be  assigned  to  different  due  dates.  Theorem  2  in  Section  5  gives  a 
computational  procedure  to  find  P^.  It  was  shown  that  large  computational  savings  can  be  realized 
by  our  procedure  compared  to  a  complete  search. 
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Appendix  ,        ,-,,    ..  • 

This  appendix  states  and  proves  several  important  results. 

Lemma  1:  Let  D,  o  and  I  be  known.  Let  I  be  such  that  there  are  three  successive  jobs,  a,  b,  and  c 
(Ca  <  Cb  <  Cc)  in  a  with  a,  c  e  Ij  and  b  2  Ij.  Then  it  is  possible  to  find  an  alternate  I  (without 
changing  D  and  a)  with  the  same  or  a  lower  cost  such  that  either  a,  b  e  Ij  and  eg  Ij  or  b,  c  g  Ij 

and  a  g  Ij.  <  .  •.   , 

Proof.  Let  be  Ik,  k  ^  j.  We  provide  a  proof  by  considering  the  following  4  cases.  Ca,  Cb  and 
Cc  are  the  completion  times  for  jobs  a,  b  and  c,  respectively.  Dj  is  the  due  date  for  jobs  in  Ij  and 
Dk  is  the  due  date  for  jobs  in  I^.  We  have  Ca  <  Cb  <  Cc. 

Case  1:  Job  b  early  in  Ik  and  early  if  assigned  to  Ij,  or  Cb  ^  min(Dj,Dk) 
A(b)    =    COST  of  bin  Ik -COST  of  b  in  Ij 

=    [PrDk+P2-(Dk-Cb)]  -  [PrDj+P2-(Dj-Cb)] 
=    Pr(Dk-Dj)+P2-(Dk-Dj) 
In  this  case,  since  Ca  <  Cb,  Job  a  is  early  in  Ij  and  will  be  early  if  assigned  to  Ik- 
A(a)    =    COST  of  a  in  Ik -COST  of  a  in  Ij 
=    Pi.(Dk-Dj)-KP2-(Dk-Dj) 
=    A(b) 
Assigning  Job  a  to  Ik  and  Job  b  to  Ij  gives  the  same  cost. 

Case  2:  Job  b  tardy  in  Ik  and  tardy  if  assigned  to  Ij,  or  Cb  ^  max(Dk,Dj) 

A(b)    =     [Pi-Dk-HP3-(Cb-Dk)]  -  [PrDj+P3-(Cb-Dj)] 
=    Pr(Dk-Dj)  +  P3-(Dj-Dk) 
In  this  case  Job  c  is  tardy  in  Ij  and  will  be  tardy  if  assigned  to  Ik- 

A(c)    =    Cost  of  C  in  Ik  -  Cost  of  C  in  Ij 
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=    [PrDk+P3-(Cc-Dk)]  -  [PrDj+P3-(Cc-Dj)] 

=    ^b 

Thus,  Job  b  can  be  assigned  to  Ij  and  Job  c  can  be  assigned  to  Ik  without  increasing  the  total  cost. 

Case  3:  Job  b  tardy  in  Ik  but  early  in  Ij,  or  Dk  <  Cb  ^  Dj 

A(b)    =    [PrDk+P3-(Cb-Dk)]  -  [PrDj+P2-(Dj-Cb)] 
=     Pi.(Dk-Dj)+Cb-(P3+P2)-(P3-Dk+P2-Dj) 
In  this  case,  Job  a  is  early  in  Ij,  and  early  or  tardy  in  Ik- 

rPr(Dk-Dj)+P2-(Dk-Dj)  if  Ca  <  Dk 

lPr(Dk-Dj)+Ca-(P3+P2)-(P3-Dk+PrDj)  if  Ca  >  Dk 

It  is  easy  to  see  that  A(a)  <  A(b).  Thus,  Job  a  can  be  assigned  to  Ik  and  Job  b  can  be  assigned  to  Ij 

without  increasing  the  total  cost. 

Case  4:  Job  b  early  in  Ik  but  tardy  in  Ij,  or  Dj  <  Cb  ^  Dk 

A(b)    =    [Pi-Dk+P2-(Dk-Cb)]  -  [PrDj+P3-(Cb-Dj)] 
=     Pi-(Dk-Dj)-Cb-(P2+P3)+(P2-Dk+P3-Dj) 
In  this  case.  Job  c  is  tardy  in  Ij  and  early  or  tardy  in  Ik. 

/(Pl-PsHDk-Dj)  ifCc>Dk 

A(c)    =    \ 

lPr(Dk-Dj)-Cc-(P2+P3)+(P2-Dk+P3-Dj)  if  Cc  <  Dk 

Since  A(c)  <  A(b),  Job  c  can  be  assigned  to  Ik  and  Job  b  can  be  assigned  to  Ij  without  increasing 

the  total  cost. 

This  completes  the  proof. 

Lemma  2:  The  TP(D,I,a)  problem  has  an  optimal  schedule  with  zero  machine-idle  time. 
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Proof.  Consider  a  schedule  in  which  the  machine  is  idle  in  the  interval  (t,t+5)  and  Job  Q  starts  at 
t+5.  We  show  that  an  alternate  feasible  schedule  can  be  constructed  such  that  there  is  zero  idle  time 
in  the  interval  (t,t+5)  and  the  total  cost  does  not  go  up. 

We  complete  the  proof  by  considering  the  following  two  cases. 

Case  1:  Job  J?  is  tardy;  i.e.,  Dp  <  t+6+tp. 

Revise  the  schedule  by  moving  Job  P  backward  by  e  =  min(5,t+5+tp-Dp).  Job  Q  starts  at 
t+5-e  in  the  revised  schedule.  Either  the  machine-idle  time  in  the  interval  (t,t+5)  goes  to  zero  or 
Job  Q  completes  on  time  (or  both)  in  the  revised  schedule.  The  cost  of  the  revised  schedule  is 
lower  because  the  tardiness  of  Job  Q  has  gone  down. 

Case  2.  Job  P  is  early  or  on-time;  i.e..  Do  >  t+b+tQ. 

Let  Job  P  be  in  position  n.  Revise  the  schedule  by  moving  the  jobs  in  positions 
n,  n-i-l,...,N  backward  by  5.  Also,  reduce  the  due  dates  of  all  these  jobs  by  6.  The  revised 
solution  costs  less  because  it  has  lower  earliness  and  due  date  penalties. 

This  completes  the  proof. 

Lemma  3:  Let  D,  I,  and  a  be  known.  Let  D  be  such  that  there  is  a  Dj  that  does  not  coincide  with 
the  completion  time  of  any  job  in  Ij  then  it  is  possible  to  find  an  alternate  Dj  without  increasing  the 
total  cost. 

Proof:  Let  Tj  be  the  number  of  tardy  jobs  and  Ej  be  the  number  of  early  jobs  in  Ij.  Let  a  and  b  be 
two  adjacent  jobs  in  a  such  that  Ca  <  Dj  <  Cb-  A  is  a  small  quantity  such  that  Dj-A>Ca  and 
Dj-i-ACb-  Let  A+  denote  the  increase  in  the  total  cost  for  jobs  in  Ij  if  Dj  is  increased  by  A  and  let  A" 
denote  the  increase  in  the  total  cost  if  Dj  is  reduced  by  A.  We  have 

A+  =  nj-A-Pi-i-Ej-A-P2-Tj-A-P3,  and 

A-  =  -nj-A-Pi-Ej-A-P2+Tj-A-P3  =  -A+. 
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If  A"*"  >  0,  we  can  decrease  Dj  upto  Ca  without  increasing  the  total  cost.  If  A-  >  0,  we  can  increase 
Dj  upto  Cb  without  increasing  the  cost.  This  completes  the  proof. 

Proof  of  Lemma  5: 

We  prove  the  lemma  for  i  =  1  and  2.  The  proof  can  be  easily  extended  for  i  =  3,  4, ...,  N. 


Proof  for  WTT(n)  >  WTiM  y^l    ' 


Consider  the  problem  with  m  =  2  and  t]  =  t2  =  ...  =  In  =  t.  (Since  WTj(n)'s  are 
independent  of  processing  times,  it  is  sufficient  to  prove  this  result  for  ti  =  t2  =  ...  =  tN  =  t.)  From 
(8),  we  get 

COST  ^     ^ 

2 


WTi(n)  =  ^^^,andWTi 


2 


To  prove  WTi(n)>WTi 


frN\^^ 


,  it  is  sufficient  to  show  that 


J 


COST(n)  >  COST|  ["y^    ' 
2 


(Al) 
for  1  <  n  <  i^T.  Note  that  COST(P)  here  denotes  the  minimum  total  cost  for  the  N-job  problem 


with  m  =  2  due  dates,  ti  =  t2  =  ...  =  tN  and  Q  jobs  assigned  to  Di.  We  provide  a  proof  by 
construction. 

Consider  a  COST(n)-solution  with  n  <  N-n.  Let  E  (T  )  denote  the  number  of  early  (tardy) 
jobs  in  Batch  1  and  E-(T  )  denote  the  number  of  early  (tardy)  jobs  in  Batch  2  in  the  COST(n)- 
solution.  It  is  easy  to  see  that  n  =  E  +T  ,  Di  =  E  -t,  D2  =  (E  +E  +T  )-t. 
Due  date  penalty  for  the  first  n  jobs  =  n-D  1  -P 1 

=  (Ej+TJ).(EJ-t)-Pi, 
Earliness/Tardiness  Cost  for  the  first  n  jobs      =  P2-((E  -l)-t+(E  -2)-t  +...+l-t) + 

P3(t+2t+...+T'-t) 
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/  / 


P2-E^.(Ej-l).t     P3-Tj.(Tj  +  l)-t 


2 
Due  date  penalty  for  the  last  N-n  jobs  =  (N-n)*D2-Pi, 


/  /  / 


=   (E^+T^)'(E^+E^+T^)'f?i. 


Pl-E^-CE^-D-t     P3-T2-(T2+l)-t 


Earliness/Tardiness  Cost  for  the  last  N-n  jobs  =  j^ + ^ ,  and 


/  / 


/  /  / 


Ej-(Ej-l)-t-P2     Tj.(Tj  +  l)-fP3 


COST(n)    =    (E^+T^HE^.t>Pi  + ^ + 2 "^ 


/  /  f  f 


/  /  f  r  f 


E2-(E2-l)-t-P2     T^-CT^+D-t.Ps 


(E2+T2)-(Ej+T^+E2)-t-Pi  + ^ + 2 " 

Now  consider  another  feasible  solution  constructed  from  the  COST(n)-solution  by 
changing  the  assignment  for  the  first  job  in  Batch  2.  We  assume  that  the  assignment  for  the  first 
job  in  Batch  2  is  changed  to  Batch  1.  Thus  the  first  (n+1)  jobs  are  now  assigned  to  Batch  1  and 
the  last  N-(n+l)  jobs  are  assigned  to  Batch  2.  The  due  dates  Di  and  D2  are  determined  as  follows. 


CASE  1:  Di  =  (E^+l)-t  and  D2  =  (Ej+E2+Tj)-t 


Let  COSTi  denote  the  cost  for  this  case,  we  have 


/  / 


(Ej  +  l)-Ej-fP2     Tj-(T^  +  l)-fP3 


COSTi    =     (Ej+Tj  +  l)-(E^+l)-fPi  + 


/  / 


/  / 


(E2-l>(E2-2).fP2     T2-(T2+l)-t-P3 


(E^+T^-\y(E^+T^+E^)'t'Pi  + 2 "^ 2 '^"^ 

COST(n)  -  COST]  =  f(Pi+P2)'(E2-Ej-l)  (A2) 


CASE  2:  Di  =  E^-t  and  D2  =  (E^+T^+E^+U-t 

Let  COST2  denote  the  cost  for  this  case,  we  have 


Ej-(Ej-l)-fP2     (T^  +  l)-(Tj-h2)-fP3 


COST2    =     (E^+Tj  +  l)-Ej-t-Pi  + 
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E2-(E2-l)-t-P2      (T^-D-T^-t-Ps 


(E^+T^-lHEj+Tj+E^+D-t-Pi  + 


,  and 


COST(n)  -  COST2  =  (T2-Tj-l)-(P3-Pi)-t 


(A3) 


Since  n  <  N-n,  we  have  either  E   <  E-  and  T  <  T^,  or  E   <  E-  and  T  <  T..  If  E   <  E  , 
then  COST(n)  >  COST],  and  if  TJ  <  T^,  then  COST(n)  >  COST  2.  Thus,  for  n  <  N-n,  an 
alternate  feasible  solution  with  an  equal  or  lower  cost  can  be  constructed  by  reassigning  a  job  from 
Batch  2  to  Batch  1.  This  implies  WTi(n)  >  WTi(n+l)  for  n  <  N-n.  This  completes  the  proof  for 


WTi(n)>WTi 

• 

Proof  for  WT?(n)  >  WTo 

lb  J  J 

Consider  the  problem  with  m  =  2,  ti  =  0  and  t2  =  13  ...  =  tN  =  t.  From  (8),  we  get 
COST(P) 


WT2(i?)  = 


•for  P  €  <1,N-1>. 


frN\^\ 


To  prove  WT2(n)  >  WT2    y      for  n  g  ( 1 ,  y    ),  it  is  sufficient  to  prove  that 

COST(n)  >  COSTJ  ry^j    |  for  n  g  <  1  ,ryY). 
Again,  a  proof  by  construction  is  provided.  Consider  a  COST(n)-solution  with  E  ,  T  ,  E.  and  T^ 
as  defined  above.  The  job  with  zero  processing  time  is  counted  in  E  .  We  have  Di  =  (E  -l)'t  and 

/  /  / 

D2  =  (E +T  -i-E--l)*t.  The  zero  processing  time  job  completes  at  Di.  We  get 


COST(n)    =    (Ej-KTjHEj-l).fPi-H 


(Ej-lHEj-2)-fP2     (Tj-^l).T^-fP3 

7^ + ^ + 


E2-(E2-l)-t-P2     T2-(T2+l)-t-P3 


(E2+T2)-Pr(Ej+T^-HE2-l)-t  + 


Now  consider  an  alternate  feasible  solution  such  that  Batch  1  has  the  first  n-i-1  jobs  and  Batch  2  has 
the  last  N-(n-i-l)  jobs  in  the  sequence  for  the  COST(n)-solution.  The  zero  processing  time  job 
completes  at  Di  where  Di  and  D2  are  determined  as  follows: 
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CASEl:  Di=Ej-tanclD2  =  (Ej+Tj+E2-l)-t.  We  get 


Ej.(Ej-l>fP2     Tj.(Tj  +  l).fP3 


COSTi     =     (Ej+Tj  +  lHEj-t)-Pi  + 2 + T 


/  / 


(E2-l)-(E2-2)-fP2     T^^VD-t-Ps 


(E^+T^-D-CEj+Tj+E^-n-t-Pi  + 2 + 2 '  ^^ 

COST(n)  -  COSTi  =  (E^-EJ-D-t-Pi  (A4) 


/  /  / 


CASE  2:  Di  =  (E^-l>tandD2  =  (Ej+T^+E2)-t.  We  get 


/  / 


(Ej-l)-(Ej-2)-fP2     (T^+l)(Tj+2)-fP3 


COST2    =     (E^+Tj  +  lHEj-l)-fPi  + ^ + ^ + 

E2-(E2-l)-t-P2      (T^-D-T^-t-Ps 


/  /  / 


(E2+T2-l)-(Ej+Tj+E2)-t-Pi  + ^ + 2 '  ^"^ 

COST(n)  -  COST2  =  (T^-Tj-D-CPs-PO-t  (A5) 

Note  that  (A4)  is  the  same  as  (A2)  and  (A5)  is  the  same  as  (A3).  As  in  the  WTi(-)  case, 
we  can  argue  WT2(n)  >  WT2(n+l)  for  n  <  N-n.  This  completes  the  proof  for  WT2(0- 
The  above  proof  can  be  extended  to  show  WTi(n)  >  WTi(n+l)  for  n  <  N-n  and 
i  e  <3,N>.  To  prove  it  for  i  =  k,  it  will  be  assumed  that  the  2-due  date  N-job  problem  has  (k-1) 
jobs  with  zero  processing  times,  and  N-(k-l)  jobs  with  processing  time  of  t  >  0  each.  All  the  zero 
processing  time  jobs  are  assumed  to  complete  at  D].  For  a  given  COST(n)  solution,  an  alternate 
feasible  solution  is  constructed  by  increasing  the  number  of  jobs  assigned  to  Batch  1  from  n  to 
n+1.  Two  alternatives  for  due  dates  are  considered.  Compared  to  the  due  dates  in  COST(n),  the 
due  date  for  Batch  1  is  increased  by  t  while  the  due  date  for  Batch  2  is  kept  the  same  in  Case  1 .  In 
Case  2,  the  due  date  for  Batch  1  is  kept  the  same  while  the  due  date  for  Batch  2  is  increased  by  t. 
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